Compactlike kinks and vortices in generalized models 
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This work deals with the presence of topological defects in k-field models, where the dynamics 
is generalized to include higher order power in the kinetic term. We investigate kinks in (1,1) 
dimensions and vortices in (2,1) dimensions, focusing on some specific features of the solutions. In 
particular, we show how the kinks and vortices change to compactlike solutions, controlled by the 
parameter used to introduce the generalized models. 
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I. INTRODUCTION 

In this work, we deal with defect solutions in k-field 
models, which are models where the kinematics is gen- 
eralized to allow for the presence of terms depending on 
higher order power of the derivative of the fields. 

There are several distinct motivations to study defect 
structures in high energy physics. In the case of standard 
models, interest to investigate defect structures can be 
found, for instance, in [HHf. Motivations to study defect 
structures in generalized models come from Cosmology, 
with the k-essence models and from other areas, as 
one can find in the recent works 

A nice property of generalized models is that under 
specific conditions, they may support compactons [l9j . 
which are defect solutions which live in a compact region, 
so they have finite wavelength. This is different from the 
standard defects, which are described by solutions of in- 
finite wavelength. Since compactons have gained recent 
interest in high energy physics [2(J, the main purpose of 
the present work is to investigate the behavior of the de- 
fect solutions, and their modification into compactons, 
under the variation of the driving parameter, which re- 
sponds for the generalized model. As we show below, the 
generalized models which we will consider are controlled 
by a single real parameter, labeled a, which responds for 
the generalization, in the sense that the limit a — > leads 
us back to the standard model. This driving parameter a 
is then used to make the generalized model close to (for 
a small) or far away from (for 1/a small) the standard 
model. 

As one knows, the existence of excitations localized un- 
der the presence of nonlinear interaction has long been 
explored with the hope to better understand the funda- 
mental contents of matter [HQ. In the case of com- 
pactons [Til [2(| , which are excitations characterized by 
having a compact support, one notes that two adjacent 
compactons do not interact unless they come into close 
contact. This is specific to compactons and in this sense, 
compact excitations seem to be well appropriate to in- 
troduce new features as particlelike structures, as kinks 



in the line or vortices in the plane, or immersed in space 
as domain walls or cosmic strings, respectively [lMl8l|. 

In the present work we deal with kinks in generalized 
models described by a single real scalar field 4> in (1, 1) 
space-time dimensions. This is done in the next Sec. HH 
where we start with the standard model and then gener- 
alize it and study the presence of defect structures. Since 
the case of kinks is simpler, we use it to set the focus of 
the work, to prepare for the study of vortices, which is 
done in Sec. lIIII We deal with vortices considering gener- 
alized models in (2, 1) space-time dimensions, described 
by a complex scalar field <p coupled to the £7(1) gauge 
field Afj,, with the standard model being the Maxwell- 
Higgs model firstly investigated in [2l|, with the gener- 
alization being controlled by the parameter a, in a way 
similar to the case of kinks considered in Sec. |TTJ We end 
the work in Sec. IIV1 where we introduce our comments 
and conclusions. 



II. THE CASE OF KINKS 

The model which describes a single real scalar field 
in (1,1) space-time dimensions is given by the Lagrange 
density 



C = Md»4>-V(<t>), 



£ = X~ V(<P), 



where we have set 



The equation of motion for static solution is 

d 2 (/> _ dV 
dx 2 d(f> 

It can be integrated to give 
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4>' 2 = V{cj>). 



(1) 
(2) 

(3) 
(4) 
(5) 
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If we choose specific potential, we can obtain topological 
solution. As a nice model, let us take the potential 



(6) 



In this case, we are using units such that both the field 
and coordinates are dimensionless. Here we obtain the 
solution 



(x) = tanh(ir), 



(7) 



after choosing its center to be at the origin x = 0. The 
energy of the solution is E — 4/3. As we know, the width 
Wk of the defect can be written in terms of the height the 
potential between the two minima. Thus, we can write 



Wk oc V 2 (0). 



(8) 



This shows that as higher the barrier between the two 
minima is, as thinner the width of the defect solution is, 
when one fixes the distance between the minima. 
Let us now consider the generalized model 



C = X - aX 2 



(9) 



where a is a real parameter, introduced to control the 
modification of the standard kinematics, with the limit 
a — > leading us back to the standard situation. Here 
the equation of motion for static solution is 



(1 - 3a, 



It can be integrated to give 



or 



¥ 2 - u - 



where we have set 



iyr+T2^)-i 

m - g — 



(10) 



(11) 



(12) 



(13) 



Note that the minima of the potential are given by 
U(4>) = 0, which imposes that V((f>) = 0, as it happens 
to be the case in the standard model. This shows that in 
the above description, the distance between the minima 
does not depend on a. However, for the potential given 
by Eq. ©, we get that 



6a - 1 



a 



(14) 



and so the width of the defect in the generalized model 
depends on a. In fact, as we have investigated numer- 
ically, the dependence of the width on a is strong, and 
may make it difficult to understand how to get to the case 
of compacton solutions. However, to make the study of 



the behavior of the defect structure for large a easier to 
understand, it is better to modify the potential V(4>) in a 
way such that U(0) is constant, independent of a. Thus, 
we consider the simple case in which U(Q) = 1/2. We get 
to this with the modification 



V{<t>)^{l + -a)V{ct>), 



(15) 



which leads to 



Wl + 6a(2 + 3q W )-l - 
6 a 

With this at hand, it is now easy to understand the 
behavior of the solution for increasing a. We do this by 
first expanding the newer £/(</>) in terms of a" 1 . In the 
case of a — ¥ oo we get to 



{/(</>) = i| W 2 1 . 
We use this into Eq. (fT2"j) to get to the compacton 



cf>{x) 



1 

sin(x) 
1 



for x > | , 



for - f 
forx < 



< x < 



(17) 



(18) 



In Fig. [T] we plot both the kink ([7]), with dashed line, 
and the compacton ([T5]). with solid line, to show how the 
kink should behave in the limit of very large a. See also 
Fig. 2 for the other plots shown in Fig. 1. 

We also depict in Fig. 3 the energy density p(x) of the 
kink- like solution for a = 0, 1, and 10. We see from this 
figure that the energy density is localized, irrespective of 
the value of alpha, as expected. 




FIG. 1: The standard kink (0) is shown with the dashed line, 
and the compacton (I18[l is shown with the solid line. The 
other solutions are for other values of a, as we explain in 
Fig. El 

In the generalized model, the first-order equation 
which we have to solve for a arbitrary is given by 



b , 2 _ 1 v/l + 3a(2 + 3a)(l-0 2 ) 2 
3 a 



1 



(19) 
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We numerically investigate this equation for several val- 
ues of a, and we also plot some solutions in Fig [T] the 
dotted line is for a = 1, and the dash-dotted line is for 
a = 100. The behavior of the solution for varying a is 
better seem in Fig. [2J There we see that for increasing a 
the defect converges to the compacton shown in Fig. ((TJ) , 
as expected. However, the convergence is very slow since 
one needs a very large a to make the kink to behave as 
a compacton. 
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FIG. 2: A closer view of Fig. 1, showing the standard kink 
(dashed line), the compact solution (solid line) and two other 
solutions, for a = 1 (dotted line) and for a = 100 (dash- 
dotted line). 

The above investigation help us to understand the ba- 
sic behavior of the defect solutions, and the passage to 
compactons. It will help us to understand the much more 
complicated situation, where we deal with vortices, to in- 
vestigate the presence of compactlike vortices in general- 
ized models. This is the subject of the next Sec. IIIII 

III. THE CASE OF VORTICES 

Let us first introduce our model. It is described by the 
Lagrange density 

£ = ™ V" + \D^\ 2 - a \D^f - V-(|0| 2 ) , (20) 

where F^ v = d lL A v - d v A^ and 

= + ieArf, (21) 

and y(|0| 2 ) is the potential which implements sponta- 
neous symmetry breaking. 

As before, here a is a real parameter which controls 
the generalized dynamics, with a — > leading us back 
to the standard Maxwell-Higgs model. Also, e stands for 
the electric charge, and the potential usually has two pa- 
rameters, A and v, which represent the coupling constant 
for self-interaction of the scalar field and the spontaneous 
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FIG. 3: The energy density p(x) of the kink-like solution of 
Eq. (|19[) for a = 0, 1 and 10, with greater a corresponding to 
higher solid curve. 



symmetry breaking parameter, respectively. Usually, it 
is given by 

V-(0) = ^V-H 2 ) 2 . (22) 

However, we are working in (2, 1) space-time dimensions, 
and so the fields, parameters and coordinates are not di- 
mensionless quantities. For simplicity, however, we can 
rescale fields, parameters and coordinates in order to 
work with dimensionless quantities. We do this in the 
usual way, and we can write the new Lagrange density in 
the simpler form 

C = -\F^ + \D^\ 2 - a \D^ - V{\<P\ 2 ) , (23) 

where we are now using dimensionless fields, coordinates 
and a, with A = v = e, for simplicity. We are also using 
the same F^ v = d^A^ — d v A^, but now we have changed 
the covariant derivative to 

D^ = d^ + iA^. (24) 

In this case, we can consider the potential in the form 

^) = J(l + ^)(l~H 2 ) 2 - (25) 

Note that the potential has the same form used in the 
case of kinks, and the presence of a has the same mo- 
tivation there considered. Note also that we could have 
introduced other models, one of them with another extra 
a-dependent contributions added to the Maxwell term, 
for instance. However, we have implemented the above 
modification since it leads to the simplest model which 
can be obtained starting from the scalar field model used 
in the former Sec. |TTJ 

In the investigation below, for simplicity we will some- 
times use 

Y = \D^\ 2 (26) 
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and 



and 



F(Y) = Y-aY 2 



(27) 



to describe the generalization included in the above 
model in a shorter way. 

The equations of motion are now given by 



and 



where 



d a F» a = , 



F Y D^D^ + FyyD^Y = 



Fy 



with F Y = dF/dY and F YY = d 2 F/dY 2 . 
The Gauss Law is written as 

V 2 A° = 2A°\(I>\ 2 F Y , 



(28) 



(29) 



(30) 



(31) 



which allows that we use A = as a proper gauge choice. 
We fix this temporal gauge from now on. 

In order to investigate the presence of vortices, let us 
consider the standard static and rotationally symmetric 
Ansatz, which implies that 



A(r 



r,6) = g (r) exp (inO) , 
1 



(a(r) — n) i 



(32) 



(33) 



where n — ±1, ±2, ... describes the vorticity of the solu- 
tion. In this case, the dimensionless static and rotation- 
ally symmetric equations of motion become 

d 2 a 1 da r ^ , . „ 9 ,„ 
d^-rTr- AaKlM = 29a > (34) 

and 

£ + + -(. + ?«)<✓-.). 

(35) 

where K±(g, a) and K2(g,a) are given by 

K 1 {g,a)=g^a(f\ + (36) 



and 



K2{g,a) 



r \ dr 



dg\ 2 d 2 g g 2 a 2 d 2 g 



dr J dr 2 



dr 2 



gadg fdg | da ga\ g"ar 



r 2 dr \ dr 

In the limit a — > 0, the above equations get to 
d 2 a 1 da 



dr 2 



dr 



2g z a 



(38) 



d 2 g ldg 
dr 2 r dr 



2 

org 



9~9- 



(39) 



which exactly reproduces the equations of motion of the 
standard Maxwell-Higgs model. According to our con- 
ventions, here we are dealing with scalar and vector fields 
with the same mass. Thus, we can write the first order 
equations 



and 



dg i 
— = ±-ga 
dr r 



(1 

r dr 



(40) 



(41) 



Their solutions are BPS states, since they solve the equa- 
tions of motion and have energy minimized to the Bogo- 
mol'nyi bound. This case is well-understood and can be 
found, for instance, in Ref. [l|. 

The generalized model is much more complicated. To 
understand the main features for a non-vanishing a, let 
us consider the energy-momentum tensor. It has the form 







= - Vltu C - F„ a F v a + 2F Y D u <t>~D^ , (42) 

where rf v = (H ) identifies the metric signature. 

Thus, the energy density is given by 



aY 2 



3 

2 a 



-a\U<f)\ 

where B is the magnetic field and 



, 3 
l + -a 



(43) 



Y 



IDc 



V (j) - iAcj) 



g 2 a 2 



The total energy of the static solution has the form 



E = / p(r) d 2 r 



(44) 



2,t / nlr\ ^- + |D0| 2 + a|D0| 4 



(45) 



and the presence of a destroys the Bogomol'nyi bound. 
We then come to the conclusion that for the generalized 
model at hand we have to study the equations of motion, 
and this requires numerical investigation. 
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FIG. 4: The numerical solution of the equation of motion for 
a(r) for a — 0, 10 and 100, shown with the solid, dashed and 
dash-dotted lines, respectively. 



FIG. 5: The numerical solution of the equation of motion for 
g(r) for a — 0, 10 and 100, shown with the solid, dashed and 
dash-dotted lines, respectively. 



To solve the equations of motion with finite energy 
solutions, the boundary conditions on g(r) and a(r) are 
given by, near the origin 



lim<? (r) — ¥ and lima(r) 

r— >0 r— >0 



and at very large distances 



lim g (r) — > 1 and lim a (r) 



n, 



0. 



(46) 



(47) 



We have investigated the equations of motion (jMJ) and 
(1351) numerically, and the results for a(r) and g{r) are 
depicted in Figs. 4 and 5 for some values of a. 

The numerical strategy was to use the pseudospectral 
method [22| that consisted in approximating the fields 
a(r) and g(r) by 



N 



a(r) ~ a N (r) = ^ a k ipk(r), 

k=0 
N 

g(r) ~ g N {r) = ^ gkXk{r), 



(48) 
(49) 



k=0 



where N is the truncation order that dictates the number 
of modes dk and gk kept in the series, with ipk(r) and 
Xk (r) being the basis functions expressed in terms of the 
Chebyshev polynomials, and reproducing the boundary 
conditions gH]) and (|47)) . 

We have used N = 15 and followed straightforwardly 
the steps for the implementation of the pseudospectral 
method in the present case resulting in the determination 
of the modes &k and gk yielding, as a consequence, the 
reconstruction of the fields a(r) and g{r). The numerical 
results are obtained for n = 1. In Figs. 4 and 5 we 
depict the fields a(r) and g(r) for several values of a. We 
note from these figures that both a(r) and g(r) go to the 



corresponding vacuum states very rapidly, for increasing 
values of the parameter a. The results show that for 
larges values of a, the vortices behave as compactlike 
solutions, becoming constant field configurations at some 
finite distance from the origin. 

An important way to study topological structures 
requires that we search for the conserved topological 
charge, which in the case of vortices is the flux of the 
magnetic field. To see this, let us introduce the topolog- 
ical current 



(50) 



which is conserved. Thus, we can write the topological 
change density 



dAy dA x 
dx dy 



(51) 



We can use the Ansatz for A given by Eq. ([3"3"|) to see 
that J° — B, and the topological charge density equals 
the magnetic field B. Thus, the topological charge has 
the form 



Q T = / d 2 rB = $i 



(52) 



and it gives the flux of the magnetic field in the plane. 

For this reason, let us then investigate the dimension- 
less magnetic field, which is given by 



B 



1 da 
r dr 



(53) 



We use the boundary conditions on a(r) to see that 

Qt = $b= 27m, (54) 

showing that the topological charge or the flux of the 
magnetic field is conserved, and it is quantized according 
to the winding number n = ±1, ±2, ... 
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In fig. 6 we plot the magnetic field for a = 0,1,5 and 
10, and there we see that it has the appropriate feature, 
showing the compactlike behavior of the vortex as a in- 
creases. This is a nice behavior, because the magnetic 
field is gauge invariant, and it is directly related to the 
topological charge of the planar vortices. We also plot in 
Fig. 7 the energy density of the vortices for a — 0, 1, 5, 
and 10. We note that the energy density does not show 
the same behavior of the magnetic field, because it is 
localized irrespective of the value of a, as expected. 




FIG. 6: The dimensionless magnetic field is plotted for a = 
0, 1, 5 and 10, with the solid, dotted, dashed and dash-dotted 
lines, respectively. 



IV. ENDING COMMENTS 

In this work we have dealt with the presence of defect 
structures in (1,1) and in (2,1) space-time dimensions. 
In (1, 1) space-time dimensions, we focused on kinks and 
the corresponding compactlike solutions in generalized 
models, controlled by a single parameter a. In this case, 
we have verified that the presence of a makes the study 
somehow complicated, and so we have changed the po- 
tential, in order to simplify the investigation. The mod- 
ification introduced very much help us to clearly see the 
behavior of the solution for increasing a, reaching com- 
pactlike features for large values of a. 



In (2, 1) space-time dimensions, we have investigated 
the presence of vortices in a generalized model, with the 
generalization being also controlled by the single real pa- 
rameter a. The case of vortices is quite different from 
the case of kinks, but here we have solved the equations 
of motion for several values of a, and we have identified 
the compactlike behavior in the vortex solution which we 
found numerically. This seems to be a new behavior of 
the vortices, and we hope that the present investigation 
will stimulate new investigations in the field, mainly on 
the main features the compactlike vortices may engen- 
der. In particular, we need more numerical investiga- 
tions to see if the compactlike behavior which we have 
found can lead to compact vortices. In this sense, an in- 
teresting issue could be to investigate if the compactlike 
vortices obey the behavior found before for standard vor- 
tices [23|]. Another issue concerns the natural extension 
of the present work to the case of generalized model in 
the presence of the Chern-Simons dynamics [24l [25[ . 




FIG. 7: Energy density of the vortices for a = 0, 1, 5 and 10, 
with the solid, dotted, dashed and dash-dotted lines, respec- 
tively. 
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